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LONG-TERM  GOALS 

Development  of  a  eonsortium  for  hybrid-eoordinate  data  assimilative  ocean  modeling,  which  will  be 
ready  in  2003  to  address  both  the  US-GODAE  (Global  Oeean  Data  Assimilation  Experiment)  principal 
objective,  i.e.,  the  depietion  of  the  three-dimensional  ocean  state  at  fine  resolution  in  near-real  time, 
and  the  elimate  modeling  objective  of  produeing  an  unbiased  estimate  of  the  state  of  the  oeean  at 
coarse  resolution  for  long-term  climate  variability  researeh. 

OBJECTIVES 

The  primary  objective  is  the  establishment  of  a  global  real-time  ocean  forecast  system  based  on  a 
hybrid-eoordinate  ocean  model  with  sophisticated  data  assimilation  techniques  that  can  be  efficiently 
executed  on  massively  parallel  eomputers. 

APPROACH 

•  The  global  model  eonfiguration  is  the  one  adopted  by  the  Eos  Alamos  National  Eaboratory 
(LANE)  in  a  comparison  between  HYCOM  and  POP  (Parallel  Oeean  Program).  Since  POP  will 
beeome  the  oceanie  component  in  the  next  generation  of  the  National  Center  for  Atmospherie 
Research  (NCAR)  Climate  System  Model  (CSM),  this  ensures  that  HYCOM's  set-up  and  foreing 
parameterizations  eonform  to  the  latest  consensus  on  climate  modeling; 

•  Coarse  resolution  model-based  reanalysis  of  archived  observational  data  will  provide  a 
comprehensive  pieture  of  the  dynamies  and  thermo dynamies  of  the  global  oeean  during  reeent 
deeades; 

•  Expertise  on  the  model's  behavior  with  an  eddy-resolving  grid  will  be  gained  by  running  the  model 
in  basin-scale  configurations  using  lateral  boundary  eonditions  provided  by  the  global  simulations. 

WORK  COMPLETED 

a)  Release  of  HYCOM  1.0 

b)  Global  and  basin-scale  simulations 
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b)  Data  assimilation  capabilities 


RESULTS 
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Figure  1:  Winter  mixed-layer  thickness  for  HYCOM  1. 0  KPP  (upper  left  panel),  MICOM  (upper 
right  panel),  HYCOM  1. 0  KTl  (lower  left  panel),  and  HYCOM  1. 0  KT2  (lower  right  panel). 
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Traditional  vertical  coordinate  choices  [z-level,  terrain- following  (sigma),  isopycnic]  are  not  by 
themselves  optimal  everywhere  in  the  ocean,  as  pointed  out  by  recent  model  comparison  exercises 
performed  in  Europe  (DYnamics  of  North  Atlantic  MOdels  -  DYNAMO)  and  in  the  U.S.  (Data 
Assimilation  and  Model  Evaluation  Experiment  -  DAMEE).  Ideally,  an  ocean  general  circulation 
model  (OGCM)  should  (a)  retain  its  water  mass  characteristics  for  centuries  (a  characteristic  of 
isopycnic  coordinates),  (b)  have  high  vertical  resolution  in  the  surface  mixed  layer  (a  characteristic  of 
z-level  coordinates)  for  proper  representation  of  thermodynamical  and  biochemical  processes,  (c) 
maintain  sufficient  vertical  resolution  in  unstratified  or  weakly-stratified  regions  of  the  ocean,  and  (d) 
have  high  vertical  resolution  in  coastal  regions  (a  characteristic  of  terrain- following  coordinates). 

The  hybrid  coordinate  adopted  in  HYCOM  (HYbrid  Coordinate  Ocean  Model)  is  one  that  is  isopycnal 
in  the  open,  stratified  ocean,  but  smoothly  reverts  to  a  terrain- following  coordinate  in  shallow  coastal 
regions,  and  to  pressure-level  coordinates  in  the  mixed  layer  and/or  unstratified  seas.  The  hybrid 
coordinate  extends  the  geographic  range  of  applicability  of  traditional  isopycnic  coordinate  circulation 
models  (the  basis  of  the  present  hybrid  code),  such  as  the  Miami  Isopycnic  Coordinate  Ocean  Model 
(MICOM)  and  the  Navy  Eayered  Ocean  Model  (NLOM),  toward  shallow  coastal  seas  and  unstratified 
parts  of  the  world  ocean. 

HYCOM  1.0  was  recently  released  and  is  the  result  of  collaborative  efforts  between  the  University  of 
Miami  (Halliwell,  Chassignet),  the  Eos  Alamos  National  Eaboratory  (Bleck),  and  the  Naval  Research 
Eaboratory  (Wallcraft).  The  capability  of  assigning  additional  coordinate  surfaces  to  the  oceanic  mixed 
layer  gives  us  the  option  of  replacing  the  slab-type  Kraus-Turner  mixed  layer  of  MICOM  by  a  more 
sophisticated  closure  scheme  in  HYCOM,  such  as  K-Profile  Parameterization  (KPP)  (see  Earge  et  al, 
1997).  Both  mixed  layer  parameterizations  are  available  in  HYCOM  1.0. 
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Figure  2:  One-year  time  series  of  interface  depths,  temperature  and  mixed  layer  depths  in  the 
Subtropical  Mode  Water  region  for  HYCOM 1. 0  KPP  (upper  left  panel),  MICOM  (upper  right 
panel),  HYCOM  1. 0  KTl  (lower  left  panel),  and  HYCOM  1. 0  KT2  (lower  right  panel). 


The  KPP  mixing  algorithm  has  several  advantages:  It  provides  mixing  throughout  the  water  column 
with  an  abrupt  but  smooth  transition  between  the  vigorous  mixing  in  the  surface  boundary  layer  and 
the  relatively  weak  diapycnal  mixing  in  the  ocean  interior.  It  works  on  a  relative  coarse  and  unevenly- 
space  vertical  grid.  It  parameterizes  the  influence  of  a  larger  number  of  physical  processes  that  other 
commonly  used  mixing  schemes.  In  the  ocean  interior,  the  contribution  of  background  internal  wave 
breaking,  shear  instability  mixing,  and  double  diffusion  (both  salt  fingering  and  diffusive  instability) 
are  parameterized.  In  the  surface  boundary  layer,  the  influences  of  wind-driven  mixing,  surface 
buoyancy  fluxes,  and  convective  instability  are  parameterized.  The  KPP  algorithm  also  parameterizes 
the  influence  of  nonlocal  mixing  of  T  and  S,  which  permits  the  development  of  countergradient  fluxes 
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Figure  3:  Ekman  spiral  representation  in  the  North  Atlantic\ 
trade  winds  region  (HYCOM 1. 0  KPP). 

During  each  baroclinic  time  step,  the  surface  fluxes  of  momentum,  heat,  and  mass  are  all  applied  to  the 
uppermost  model  layer  with  the  exception  of  the  shortwave  radiation  flux,  which  penetrates  beneath 
the  upper  layer  with  an  exponential  decay  scale  dependent  on  the  Jerlov  water  type  (assumed  equal  to 
3).  The  KPP  algorithm  then  calculates  vertical  profiles  of  viscosity  along  with  T  and  S  diffusivity 
coefficients  at  each  interface.  These  coefficients  are  parameterized  as  functions  of  boundary  layer 
thickness  and  a  third-order  polynomial  "shape"  function,  the  latter  insuring  an  abrupt  but  smooth 
transition  between  the  small  diffusivities  of  the  ocean  interior  and  the  large  diffusivities  in  the  surface 
boundary  layer.  Given  the  profiles  of  viscosity  and  T,  S  diffusivities,  the  vertical  diffusion  equation  is 
solved  for  each  variable  using  a  tri-diagonal  matrix  algorithm.  The  solution  is  semi-implicit  and 
performed  in  two  iterations.  For  the  second  iteration,  the  vertical  viscosity  and  diffusivity  profiles  are 
re-calculated  from  the  momentum,  T,  and  S  profiles  that  were  calculated  during  the  first  iteration. 

Since  HYCOM  equations  are  solved  on  the  C  grid,  the  two  iterations  are  first  performed  on  pressure 
grid  points  using  velocity  components  interpolated  to  these  points.  The  viscosity  profiles  calculated  on 
the  pressure  points  are  then  interpolated  to  the  u  and  v  grid  points  where  the  vertical  diffusion  equation 
is  solved  once.  Mixed  layer  thickness  is  diagnosed  as  the  interpolated  depth  at  which  the  density  first 
exceeds  the  surface  layer  density  by  a  specified  amount.  Mixed  layer  variables  are  calculated  by 
integrating  them  over  the  depth  of  the  mixed  layer. 
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The  Kraus-Turner  (KT)  mixed  layer  model  was  also  implemented  in  HYCOM  1.0  as  a  benchmark  for 
comparison  with  MICOM  and  for  comparison  with  new  mixed  layer  models  such  as  the  KPP.  The 
implementation  of  the  KT  mixed  layer  model  differs  substantially  from  MICOM  since,  in  HYCOM, 
the  prognostic  mixed  layer  thickness  does  not  coincide  with  any  vertical  coordinates  and  the  mixed 
layer  base  cuts  across  layer  interfaces.  To  calculate  the  evolution  of  mixed  layer  thickness  in  the  KT 
model  and  entrainment  fluxes  across  the  mixed  layer  base,  it  is  necessary  to  quantitatively  estimate  the 
jumps  in  momentum  and  thermodynamical  variables  that  exist  at  the  mixed  layer  base.  This  is  trivial  in 
MICOM  since  the  interface  at  the  bottom  of  layer  one  is  the  mixed  layer  base.  For  HYCOM,  an 
"unmixing"  algorithm  was  developed  to  estimate  these  jumps  at  each  grid  point  by  dividing  the  layer 
containing  the  mixed  layer  base  into  two  sublayers.  The  upper  sublayer  is  part  of  the  mixed  layer  while 
the  lower  sublayer  is  part  of  the  stratified  interior.  Since  the  KT  algorithm  governs  mixing  only  within 
the  mixed  layer,  a  separate  diapycnal  mixing  algorithm  must  be  included  for  the  interior  ocean. 
HYCOM  1.0  has  been  equipped  with  two  diapycnal  mixing  algorithm  choices:  the  one  traditionally 
used  in  MICOM  (KTI)  and  an  implicit  procedure  (KT2)  that  is  essentially  the  interior  diapycnal 
mixing  component  of  the  KPP  mixing  algorithm. 

The  performance  of  HYCOM  1.0  were  compared  to  MICOM  2.8  in  a  2°  North  Atlantic  configuration 
with  monthly  climatological  COADS  forcing.  Three  primary  cases  were  run  using  HYCOM  1.0:  KPP 
mixing,  KT  mixing  with  the  explicit  diapycnal  mixing  used  in  MICOM  2.8  (KTI),  and  KT  mixing 
with  the  implicit  diapycnal  mixing  extracted  from  the  KPP  mixing  algorithm  (KT2).  The  latter 
diapycnal  mixing  algorithm  parameterizes  the  influence  of  background  internal  wave  breaking,  double 
diffusion,  and  shear  instability.  The  model  is  first  spun  up  for  thirty  years  using  the  COADS 
climatological  annual  cycle  forcing,  then  run  for  one  year  with  three-dimensional  fields  saved  monthly 
for  further  analysis.  Time  series  of  all  variables  were  also  saved  at  high  time  resolution  (five  times  per 
day)  at  locations  representing  dynamically  distinct  regimes. 

Model  performance  for  the  four  cases  is  illustrated  by  the  maps  of  winter  (15  March)  mixed  layer 
thickness  (Fig.  1)  and  by  1-year  time  series  of  mixed  layer  thickness,  interface  depths,  and  temperature 
in  the  upper  800  m  at  a  grid  point  in  the  Subtropical  Mode  Water  (STMW)  formation  region  (Fig.  2). 
The  winter  mixed  layer  depth  (Fig.  1)  shows  qualitatively  similar  patterns  for  the  four  cases.  The 
mixed  layer  over  most  of  the  subtropical  gyre  is  less  deep  in  HYCOM  1.0  KPP  than  in  MICOM  2.8 
and  in  the  two  HYCOM  I.O  cases  using  the  KT  mixed  layer  model.  The  KPP  mixing  model  tends  to 
produce  a  thinner  winter  mixed  layer  in  this  region  than  the  KT  model.  The  mixed  layer  is  also 
shallower  in  the  region  west  of  central  Europe  in  the  two  HYCOM  1.0  KT  cases  than  in  the  other  two 
cases.  The  annual  cycle  of  temperature  and  mixed  layer  thickness  in  the  STMW  formation  region  (Fig. 
2)  also  displays  qualitatively  similar  behavior  among  the  four  cases  with  the  deepest  winter  mixed 
layer  found  in  the  MICOM  simulation.  High  vertical  resolution  in  the  HYCOM  surface  mixed  layer 
now  allows  for  a  representation  of  the  Ekman  dynamics  as  illustrated  in  Eig.  3. 

Eor  assimilation  of  surface  altimetry  data,  the  Reduced  Order  Information  filter  (ROIE)  has  been 
ported  from  MICOM  to  the  HYCOM  coordinate  system.  ROIE  is  a  reduced  order  implementation  of 
Kalman  filter  that  has  been  shown  (Chin  et  al.  1999;  2000)  effective  in  observation  system  simulation 
experiments  (twin  experiments)  for  assimilation  of  sea  surface  heights  sampled  sparsely  as  in  the 
Topex/Poseidon  satellite.  The  assimilation  technique  has  also  been  proven  mathematically  to  preserve 
positive  definiteness  in  the  approximated  covariance  matrix  (Chin,  2000). 

Porting  of  ROIE  assimilation  system  from  MICOM  to  HYCOM  is  achieved  by  converting  its  vertical 
axis  from  one  that  based  on  layer-thickness  to  a  pressure-based  coordinate.  The  HYCOM-adapted 
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ROIF  is  tested  first  in  a  2  x  2  degree  grid  model  of  North  Atlantic  with  the  KPP  mixing  scheme. 
Preliminary  analysis  show  that  the  sea  surface  height  field  converges  exponentially  towards  the  truth  (a 
twin  run)  through  assimilation  of  the  sparsely  sampled  (Topex/Poseidon  track)  data.  This  result  is 
consistent  with  the  previous  performance  with  MICOM.  The  data  were  assimilated  at  every  model 
time  step  (if  available),  instead  of  every  10  days. 

Further  examination  of  the  HYCOM-adapted  ROIF  is  underway.  Our  agenda  is  to  perform  (i)  quality 
checks  of  the  new  ROIF  numerics,  e.g.,  matrix  conditioning,  (ii)  an  evaluation  of  the  vertical  profiling 
and  its  sensitivity  to  the  prior/empirical  statistics,  (hi)  an  evaluation  with  finer  resolution  HYCOM, 
i.e.,  1/3-degree  North  Atlantic. 

A  parallel  development  is  an  examination  of  the  use  of  drifter  data  to  reduce  geoid  uncertainty 
associated  with  the  Topex/Poseidon  altimetry  data.  Initial  runs  of  ROIF-based  assimilation  of  drifter- 
based  velocity  data  have  recently  completed.  The  next  step  is  co-assimilation  of  these  velocity  data 
and  altimetry  data. 

IMPACT/APPLICATIONS 

Generation  of  optimal  estimates  of  the  time- varying  ocean  state  in  support  of  the  NAVY's  needs  on 
synoptic  time  scales  on  the  order  of  weeks  to  months  and  on  spatial  scales  typically  on  the  order  of  10- 
1000  km  (mesoscale). 

TRANSITIONS 

None 

RELATED  PROJECTS 

This  effort  is  part  of  a  multi-institutional  NOPP  project  which  includes  E.  Chassignet  (Coordinator),  T. 
Chin,  G.  Halliwell,  and  A.  Mariano  (U.  of  Miami/RSMAS),  R.  Bleck  (LANE),  H.  Hurlburt,  A. 
Wallcraft,  P.  Hogan,  R.  RHodes,  C.  Barron,  and  G.  Jacobs  (NRL-Stennis),  O.M.  Smedstad  and  J.E. 
Cayula  (Planning  Systems,  Inc.),  W.C.  Thacker  (NOAA/AOML),  and  M.  O'Keefe  (U.  Minnesota). 
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